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Sterile Neutrino Sensitivity with Wrong-Sign 
Muon Appearance at z/STORM 
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Figure 1: A diagram of the proposed ^STORM facility. 



Neutrinos from STORed Muonfl (VSTORM) is a proposed experiment that 
uses 3.8 GeV/c muon decay to produce a well-understood beam of electron 
and muon neutrinos that can be used for short baseline physics (Fig. [IJ. A 
magnetized far detector allows for the wrong-sign muon appearance physics of 
v e — > and provides more sensitivity to sterile neutrinos than other proposals 
(See comparisons in [1]) because of the relative ease to which muon tracks can 
be identified. Other physics such as v e and cross section measurements are 
possible. For further details, see the Letter of Intent [16]. 

An explanation of the j/STORM appearance analysis will follow. This work 
is a continuation of the work presented in |20) . For disappearance measurement 
work, see (21"] . 



1 Short Baseline Oscillations 

LEP experiments revealed that there are three light neutrinos that couple to the 
Z-boson (ie. active neutrinos) , however, there are theoretical and experimental 
motivations [T] for neutrinos without Standard Model interactions called sterile 
neutrinos. The (3+1) scenario is the case of three active neutrinos with an 
additional heavy sterile neutrino - >> mothers ~~ and only this situation is 
considered although the results are generalizable. 

The probability v e — > depends on the mixing matrix U (Reviewed in [5]). 
Let Rij be a rotation between the i-th and j-th mass eigenstates without a CP 
violating phase: CP violation cannot be observed in oscillations with large Am 2 
dominance (See p.g. 273 of [7]). For N neutrinos, has dimension N x N 
and takes the form: 

x The facility was previously called the Very Low Energy Neutrino Factory (VLENF). 
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Table 1: Matrix elements for muon decay 
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By convention, the three neutrino mixing matrix is C/pmns = -^23^13-^12- In 
the (3+1) model of neutrino oscillations, extra rotations can be introduced such 
that the mixing matrix is Urs+i) = Rs&RzaRia ^pmns- Given that Am^ >> 
Am^, J/pmns can be approximated by the identity matrix (ie. the "short 
baseline approximation") implying U e 4 = sin(#i4) and = sin(#24) cos(#i4). 

The oscillation probabilities for appearance and disappearance, respectively, 

arc: 
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in this short baseline limit where the definition sin (28 efJ/ ) 
been introduced. 

Electron and muon neutrino disappearance measurements will constrain 
| t/ e 4 1 2 ([H]) and |f7 M 4| 2 while the appearance channel analysis could measure 
the product |C/ e 4| 2 |L r /J 4| 2 . Information about the matrix element U e 4 also arises 
from jointly analyzing za, disappearance and appearance. The remaining ma- 
trix element U T 4 can be extracted by analyzing NC rates |/7 S 4| 2 = ^ e T \U a 4\ 2 , 
using the other channels to constrain |C/ e 4| 2 and |J7 M 4| 2 , and assuming unitarity. 

2 The Neutrino Flux: $ 



Muon-decay beams contrast pion-decay beams because the beam characteristics 
and production mechanisms are well-known. The neutrino flux arises from the 
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Figure 2: The flux of v e and for a 3.8GeV/c muon-decay without oscillations 
at 2000 meters. No smearing due to accelerator effects has been performed. 



Electroweak decay of \i — > v^v^e and it is sufficient to compute matrix elements 
at tree level. The neutrino spectrum for a ^ — > e + is e (v e ) + v^iy^) decay in 
the rest frame of the muon follows: 



(5) 



where x = 2E^ om ' jm^ £ [0, 1] is the scaled neutrino energy in the rest frame, fi 
is the solid angle in the rest frame, fo(x) and fi(x) are muon decay parameters, 
and V is the polarization. Electron and neutrino masses are negligible for this 
process and ignored, hence the inclusive range for values of x. These muon decay 
paramters can be computed to leading order with Electroweak theory (See, for 
example, chapter 6 of Ref. [TB]) and are neutrino flavor dependent (See Tabled]). 

The polarization V is set to zero, similar to other studies, and has been 
shown to average to zero due to Thomas Precession. Boosting the neutrino 
distributions into the lab frame leads to: 



dydA 

d 2 N e 
dydA 



4n 

—J^E^y 2 {\ - cos 4>) [3ml - 4E»y(l - /3 cos <f>)] (6) 
(1 - /3cos0) [ml - 2E? iy (l - /3cos0)] (7) 

ttij m .. 



where y = E v j ' is the scaled neutrino energy in the lab frame, (3 — y 1 — rn/^/E^, 
A is an area, and is the number of muons. These neutrino distributions 
(Fig. [2J are for a point source so they are not directly applicable to the decay 
straight of t^STORM. 
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Figure 3: The unoscillatcd flux of v e and for a (3.8±0.38) GcV/c muon-decay 
at 2000 meters. Accelerator effects arc included; sec the text for details. 



The number of muons assumed is 1.8 x 10 18 and is based on 10 21 protons on 
target (POT) at 60 GcV/c. It corresponds to roughly 5 years of running with a 
100 kW target station. The number of useful muon decays is motivated in [T6] . 

When computing the flux for ^STORM, the far detector approximation of 
a point-source accelerator and detector no longer is applicable since the size 
of the detector and accelerator straight (150 meters) are comparable to the 
baseline of 2000 meters. The neutrino fluxes are computed by integrating over 
the decay straight, transverse beam phase space, and detector volume. The 
beam occupies a 6D phase space (x, y, z, p x , p y , p z ) and the detector has a 
5 m x 5 m cross section with the depth set by the desired fiducial mass of 1.3 
kt. Both transverse 2D phase spaces are represented by the Twiss parameters 
a = and (3 = 40 m where the lcr Gaussian geometric emittance is assumed 
to be 2.1 mm. The spread in, for example, x is o~ x = y/fie and the angular 
divergence in x is o~ x i = \J tj 'f3. The longitudinal phase space (z and p z ) is 
described by assuming a uniform distribution in z € [0, 150 m] - accurate to 
0.5% - and p z e [3.8 ± 0.38 GeV/c]. 

The flux is computed by Monte Carlo (MC) integration: random points are 
chosen within the beam phase space and within the detector volume to de- 
termine the expected flux. This integration introduces a new computational 
requirement: the baseline is a variable that affects both the oscillation probabil- 
ity (L/E) and the flux (L~ 2 geometric factor). The GLoBES software (version 
3.1.10) O QJ] that is used for neutrino factory phenomenology treats these 
as separable problems and was modified to compute this flux (and later the 
event rates and sensitivities). Specifically, GLoBES is modified such that both 
the flux and oscillation probability are computed in the oscillation probability 
engine. The code for the analysis is available [19] under the GPL license [5J. 

The resulting flux after the integration (Fig. [3]) is corrected for accelerator 
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Figure 4: The flux at the far detector for a (3.8 ± 0.38) GeV/c muon for initial 
v e states including integration over the beam envelope and detector volume. 
Final states include v e without oscillations and both v e and with best fit 
short baseline oscillations. The normalization is 10 21 POT. 

effects. The corrections are small for far detector physics (Compare to Fig. [2} 
but are important for near detector physics where the baseline is smaller than 
the decay straight . 

3 The Oscillation Probability: (Prob.) 

This section will discuss how sterile oscillation phenomenology relates to con- 
ducting the proposed experiment. For instance, for a point-source baseline of 
2000 meters, it is possible to determine the oscillation probability (Fig. [5]) using 
Eq. [2] for any combination of L and E. 

The best fit parameters for the "short baseline anomaly" and 3x3 mixing 
(i.e. sin 2 (26'i3), Ato 2 2 , etc.) are used throughout the analysis. The best fit 
parameters for the LSND anomaly come from [5] (See TABLE [2]) and agree 
with those published by the LSND collaboration [2J. For completeness, oscilla- 
tions between known mass eigenstates are included despite not influencing the 
sensitivity: the correction is order 10 -5 . The best fit data from [10] is used to 

Table 2: Best- fit oscillation parameters for the (3+1) sterile neutrino scenario 
using combined MB v and LSND v data [8]. 
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Figure 5: The oscillation probability for the "golden channel" v e — > from 
Eq. fusing the (3+1) oscillation parameters in TABLE [2] A baseline of 2000 
meters is assumed. 



specify standard 3x3 oscillations. Without loss of generality, normal hierarchy 
is assumed and the values of known 3x3 mixing can be seen in Table |3l Errors 
associated with these quantities are ignored. 

Computationally, the SNU (version 1.1) add-on [13l [14] has been used to 
extend computations in GLoBES to 4 x 4 mixing matrices. 

4 Cross section: a 

Cross sections are required for each neutrino flavor (y^, 9^, v e , u e ) and each 
interaction type (CC or NC). The nucleon cross sections (Fig. [6]) are calculated 
in [15] and [17] for the low energy and high energies, respectively. NC cross 
sections are flavor independent. The CC cross sections are approximately flavor 
independent: Fermi's Second Golden Rule results in the same matrix elements 
and, at these energies, the phase spaces for the final-state electrons and muons 
are equal. 

The total cross section requires knowing the number of nucleons in addition 
to the nucleon cross section. The fiducial mass of 1.3 kt determines the number 
of nucleons via Avogadro's number. 



Table 3: Values for 3 x 3 oscillations used. 



sin* 0i2 = 0.319 
sin 2 6» 23 = 0.462 
sin 2 0i3 = 0.010 



Ara 2 . 



|j = 7.59 x 10- 5 eV 2 
|i = 2.46 x 10~ 3 eV 2 
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Figure 6: Neutrino cross sections per nucleon. 



5 Interaction rates: A^i nt . 

The number of neutrino interactions is computed which does not require as- 
sumptions about the detector. The interaction rates can be computed by 
Ni n t. = $ x (Prob.) x a, for flux $, oscillation probability (Prob.), and cross 
section a, where all of these quantities have been computed in the previous 
sections. 

Using the LSND anomaly best fit (TABLE [2]) as a example for a sterile 
neutrino signal, the event rates for fi + and decays are shown in TABLE |H 
Various deductions can be made about these event rates and their statistical sig- 
nificance. With either stored /i + s or stored /i~s, the statistical significance of all 
channels is greater than lOcr. Combining the NC channels together results in a 
statistical significance of 20a and 17a for stored n + and p~ , respectively. There 
are no known physics backgrounds to neither v e — > v m u CC nor riu e — > v m u CC 
interactions except to negligible solar-term oscillations, so the backgrounds will 
arise from how well the detector can differentiate these interactions. 

The number of events can also be determined as a function of energy since 
the evolution of p, a, and (Prob.) as a function of energy is known. These 
distributions are shown in Fig. [7] 

There are numerous channels with reach into the sterile neutrino parameter 
space. Most other experiments have one channel to explore (See pQ for list of 
experiments), whereas in the best case j/STORM allows for 10 signals and in 
the worst case 6 (i.e. combine v e v e CC and all NC channels). 



6 Event rates after cuts 

It must be determined how many of the raw events pass analysis cuts. Similar 
analyses have been performed for Neutrino Factories exploring CP violation 
at energies ranging from 25 GeV [4] to 5 GeV [6], but never at 3.8 GeV/c. 
Preexisting experience and knowledge exists as to fractional background levels 
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Figure 7: True channel rate energy distributions assuming the LSND anomaly 
best fit values. The transitions v e — >• v^, u e D^, — > v^, and — > are 
shown. 



Table 4: Truth event rates for 10 21 POT for the no oscillations and short baseline 
oscillations described by TABLED The statistical significances are computed. 
The combined statistical significance of NC events are 20 and 17 for stored 
/i + and respectively. There are no physics backgrounds to v e — > v m u CC 
interactions. 
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and analysis difficulties; work had to be performed in order to tune the analysis 
for this energy range. 

The detector performance can be represented by migration matrices (also 
known as response matrices or energy smearing matrices) that describe both 
the energy resolution and detection efficiency. If events are binned in terms 
of true neutrino energy then the migration matrix is needed to transform the 
distribution into the space of measured neutrino energies. For example, take 
the histogram: 

jLtrue / Artrue yirtruc /irtrue \T /q\ 

— \ ly 0.0 - 0.1 GeVi iv 0.1 - 0.2 GeV) • • • J iv 3.9 - 4.0 GeV J i l°7 

where NqqZ 1 GcV is the number of events in the bin with ranges 0.0 and 0.1 
GeV. The migration matrix M used for this analysis is a square matrix and 
defined such that £ moasurod = Mh tmc where h moasurod is the expected histogram 
of reconstructed quantities in the detector. 

With a perfect detector M = diag.(l, 1, . . . , 1). M is unitary if and only if it 
describes only energy smearing. Efficiencies are included into M by removing 
the unitarity constraint. 

Migration matrices have been computed for CC, P M CC, NC, and v e CC 
(See [16]) and can be seen in Fig. [HI The background level of v e NC events into 
the signal window are negligible compared with NC due to the lower energies. 
These numbers are derived using a GENIE and Geant4 simulation, described in 
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Figure 8: Migration matrices for CC, v NC, CC, and v e CC. 
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Figure 9: The rule rate as a function of observed energy for the appearance 
channel v e — s- v^. Migration matrices are used for CC, CC, NC, and 
v e CC. 



the cited text, which are MC method softwares. Statistical fluctuations exist in 
the migration matrices due to computational limitations. 

7 Statistics 

It is necessary to determine if the number of events observed after cuts (i.e. rule 
rates) is statistically significant. The experiment must reject the null hypothesis 
when accounting for statistical fluctuations. 

The hypothesis Ho of no oscillations is the null hypothesis and designate Hi 
to be the alternate hypothesis. These hypotheses have oscillation parameters as- 
sociated with them: let 8q ~ {Am^, #34, #24, #14} be the oscillation parameters 
associated with Hq, and similarly B\ for H\. 

The test statistic X is a function of the experimental observations and let 
W be the space of all possible values of X. One can divide W into two regions: 
the region w for those possible values of X which would suggest that the null 
hypothesis Hq is not true and the remaining region W — w. 

It is desirable to have a small probability of X - by statistical fluctuations 
alone - taking a value in w when Hq is true. A level of significance a can be 
defined: 

P(X e w\H Q ) = a (9) 

where a corresponds to, colloquially, "5cr" when a - 2.8 x 1CT 7 and "IOct" 
when a ~ 7.6 x 10~ 24 . The number of "cr" correspond to the p- value of having a 
greater than no upward fluctuation of a Gaussian centered at zero. No Gaussian 
assumptions are made in this analysis. 
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The test statistic that will be used for hypothesis testing is the likelihood 
ratio test. Let there be N observations X = {Xi, Xn} and a probability 
distribution function f{Xi\0). The likelihood function is: 



/v 



L(X|0) = XYfiXW) (10) 

i=l 

= II- a a, x a,: (ii) 

i 

where A, is the expected number of background in the bin with Xi events and 
is a function of 9. The distribution is Poisson because the background levels are 
small. The short baseline parameters 9\ for H\ are free to take any value but the 
parameters 9q are fixed to zero by the null hypothesis requiring no oscillations. 
The likelihood ratio test defines a test statistic A such that: 

A= L ( X l g o) (12) 
max 9l L(X|0i) V ' 

where the denominator is maximized with respect to 6\ while the numerator 
remains fixed. Using Eq. Illl leads to: 

\ = l[e- x * +x *(\ i /X i ) Xi . (13) 

i 

The x 2 can be defined as \ 2 = —2 In A (See [3]) which is preferable to using 
A because of specifics about how multiplication is performed by a computer. 
Using this definition, one finds: 



X 2 = -21nA = 2^A i -X i + X i ln^J (14) 

which has two degrees of freedom since the numerator of Eq. [T2] has no degrees 
of freedom and the denominator has two degrees of freedom. 

A test statistic has been defined that allows for determining if an experi- 
ment is sensitive to various oscillation parameters. The \ 2 can be computed in 
terms of energy bins, with the appropriate definition of Xi, allowing for spectral 
information to be used when computing sensitivities. 



8 The Appearance Analysis 

The parameters to be explored in the appearance analysis are Am|j and s'm 2 (9 efJ ,) 
Contours in the neutrino parameter space Am^ versus sm 2 (9 etI ) can be used 
to compare the sensitivities of various proposed short baseline experiments. A 
statistics-only x 2 using spectral information is used (FigfTU)). 
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Figure 10: Sterile sensitivity under the appearance channel v e — > v^. This 
channel is the CPT of the LSND anomaly v^—l v e . There is 10er sensitivity to 
the LSND and MiniBooNc 99% confidence interval [8]. 

Care must be taken when defining ^(Amfj, sin 2 (6* eAt )) to ensure that it 
is well-defined. In the (3+1) scenario, the signal v e — > v ^ depends on the 
amplitude sin 2 (0 e/i ) = 4|£/e4| 2 |£7 M 4| 2 and frequency Am^ (See Eq. [2]). If there 
is an appearance signal, then | J7 e4 1 2 1 [7^4 1 2 7^ which implies that both £/ e4 and 
11^4 are nonzero. There is disappearance of the CC and NC backgrounds (See 
Eq.@| which affects the background estimation in the \ 2 ■ This issue is addressed 
by not oscillating the backgrounds thus overestimating the backgrounds. 

As the cuts-based detector performance section improves and various cost 
optimizations are done, there are numerous parameters that can be tuned to 
compensate and conserve the physics that can be done with such a facility. For 
example, the optimization of baseline and energy (Fig. II 1 [) allows one to change 
the baseline depending on site constraints or modify the energy of the ring 
if the accelerator gets too expensive. As the cuts-based detector performance 
improves, the various background rejections (Fig. [12] and I13[) may allow for 
a smaller detector or cheaper target station. The tools have been developed 
that allow the important accelerator and detector performance metrics into cost 
optimizations. 

Figure [TT] shows that, for a fixed baseline, increasing the muon energy is 
always advantageous. This effect arises because the maximum of the v e flux is 
not at the oscillation maximum but rather at a higher energy. At high energies 
the oscillation probability is: 
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Figure 11: A baseline optimization using a total rates statistics-only x 2 , a 
signal efficiency of 0.5, and background rejection of charge misidcntification and 
NCs at 10" 3 and 10" 4 . 

The oscillation probability decreases as E~ 2 for a fixed baseline. The signal 
rates increase as E 3 : there is a factor of E 2 from the solid angle arising from 
the I/7 opening angle and another factor of E from the cross section. The 
conclusion is that raising the stored muon energy will increase the event rates 
linearly with energy for a fixed baseline. This result has been confirmed by 
similar analyses for other muon-decay based facilities (See sensitivity work in 
i)- 

9 Conclusion 

The sensitivity of ivSTORM rules out the LSND 99% confidence interval at 
lOcr using only appearance information. The appearance channel is the CPT 
invariant of the observed anti-neutrino LSND anomaly. Optimizations have 
been shown to guide future costing and performance work. 
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Signal Efficiency 

Figure 12: Tuning the NC rejection cut. The NC rejection level is shown versus 
the signal efficiency. A charge misidcntification background of 10~ 4 is shown to 
illustrate when NC backgrounds become statistically significant. A total rates 
statistics-only \ 2 is used. 




Signal Efficiency 

Figure 13: Tuning the charge misidcntification cut. The charge misidentifica- 
tion level is shown versus the signal efficiency. A NC background of 10 -4 is 
shown to illustrate when charge misidcntification backgrounds become statisti- 
cally significant. A total rates statistics-only x 2 is used. 
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Figure 14: An optimization between the detector performance and accelerator 
performance using the charge misidentification rates and number of muon decays 
as the performance metric. IDR refers to the Interim Design Report [4j detector 
performance. FODO refers to the FODO lattice design that gives 1.8 x 10 18 
useful muon decays whilst FFAG refers to the FFAG design that gives 4.68 x 10 18 
useful muon decays. Both accelerators assume a front-end of the main injector 
at 60 GeV/c. 
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